Clustering of SCUBA galaxies and implications for the Planck mission 
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Abstract. The clustering properties of SCUBA-selected galaxies are investigated within the 
framework of a unifying scheme relating the formation of QSOs and spheroids. The theoretical 
angular correlation function is derived for different bias functions, corresponding to different 
values of the ratio -/Vfh a i /Af S ph between the mass of the dark halos hosting such galaxies and 
the mass in stars produced at the end of the major star-formation burst. SCUBA sources are 
predicted to be strongly clustered, with a clustering strength increasing with mass. Comparisons 
with the best available measurements show better fits for Mh a i /M S ph — 100. The model can 
also account for the clustering of Lyman-break galaxies, seen as the optical counterpart of low- 
to intermediate-mass primeval spheroidal galaxies. Best agreement is once again obtained for 
high values of the Mh a i /M sp h ratio. We also discuss implications for small scale fluctuations 
observed at different wavelengths by forthcoming experiments such as the Planck mission aimed 
t*"- ' at mapping the Cosmic Microwave Background (CMB). 
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1 Clustering of Scuba Sources 



The theoretical expression for the angular two-point correlation function w{0) can be derived from its 
spatial counterpart £ by projection via the relativistic Limber equation (Peebles 1980): 

O ■ 

U : Jo 00 Jo 00 N 2 (z) b 2 cS (M min , z) (dz/dx) g(r, z) dz du 
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where x is the comoving radial coordinate, r = (u 2 + x 2 6 2 ) 1 ' 2 (for a flat universe and in the small 
angle approximation), and N(z) is the number of objects within the shell (z, z + dz). 
The mass-mass correlation function £(r, z) to be inserted in eq.(|l|) has been obtained following the 
work by [10 1 (see also Q and ||]), which provides an analytical way to derive the trend of £{r,z) 
both in the linear and non-linear regime. Note that £(r, z) only depends on the underlying cosmology, 
which we fix by adopting h§ = 0.7, 0,q = 0.3, A = 0.7 and a COBE- normalized value of dg = 1. The 
relevant properties of SCUBA galaxies are included in the redshift distribution of sources N(z), and 
in the bias factor 6 e g (M m i n , z). 

The effective bias factor b e g(M m i a , z) of all the dark matter haloes with masses greater than some 
threshold M m [ n is then obtained by integrating the quantity b(M, z) (whose expression has been taken 
from P]) - representing the bias of individual haloes of mass M - opportunely weighted by the number 
density nscuBA(M,z) of SCUBA sources: 

, M _ J^ min dMb{M,z)n S cvBA(M,z) 
CS{Z) ~ IS mi JMn SC VBA(M,z) ■ [2) 

Note that, as nscuBA can be thought as the fraction of haloes hosting a galaxy in the process of forming 
stars, its expression can also be written as nscuBA(M, z) = n(M, z) Ts/th, where n(M, z) dM is the 
mass spectrum of haloes with masses between M and M + dM ( p4| ), Tg is the duration of the 
star-formation burst and th is the life-time of the haloes in which these objects reside (see ||). 

According to Q, sources showing up in the SCUBA counts can be broadly divided into three 
categories: low-mass (masses in the range M sp h — 10 9 — 1O 1O M0, duration of the star formation burst 
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Figure 1: Predictions for the angular correlation function w(9) of forming spheroids at 850 /im for 
different halo-to-bulge mass ratios and different flux cuts. Solid lines are for sources brighter than 
50 mJy, dashed lines for sources brighter than 10 mJy, while dotted lines correspond to S > 1 mJy. 
Higher curves of each type correspond to Mh a i /M sp h = 100, lower ones to Mh a i / M sp ^ = 10. 
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Tb ~ 2 Gyr, and typical fluxes S <^ 1 mJy), intermediate-mass (M sp h — 10 10 -10 11 Mq and Tb ~ 1 Gyr) 
and high-mass (M sp h ^ 10 n MQ, Tb ~ 0.5 Gyr, dominating the counts at fluxes S ^ 5 — 10 mJy). 
Note that by M sp h we denote the mass in stars at completion of the star formation process. 
In order to evaluate the bias factor in eq.(|2J) we then consider two extreme cases for the ratio between 
the mass in stars and the mass of the host dark halo: -Mhaio/^sph = 100 and Mhaio/Msph = 10. 
Mhaio/Mgph = 10 roughly corresponds to the ratio ^o/^baryon between total and baryon density, 
where we adopted for the latter quantity the standard value from primordial nucleosynthesis; this 
corresponds to having assumed all the baryons to be locked into stars and, as a consequence, has to 
be considered as a conservative lower limit. Mhaio/-^sph = 100 is instead related to S7o/T2*, f2* being 
the present mass density in visible stars: the likely value is expected to be M sp h/Mh a io = 1 — 3 %. 

Armed with the above results we can then evaluate the two-point correlation function in eq 
for different -Mhaio/Msph ratios and different flux cuts. Figure 1 presents our predictions for w(9 
respectively for a flux cut of 50 (solid line), 10 (dashed line) and 1 (dotted line) mJy. Higher curves 
of each kind correspond to the case .Mhaio / M sp h = 100, while lower curves refer to -/Vfh a io/-^sph = 10. 
The highest clustering amplitude is found for the brightest sources (S > 50 mJy). This is because 
they are associated to the most massive dark halos and are therefore highly biased tracers of the dark 
matter distribution. In addition, according to j3j, they have a rather narrow redshift distribution so 
that the dilution of the clustering signal is minimum. The very sharp drop of all the curves at 8 ~ 1° 
is due to the absence of nearby objects. This reflects the notion that the actively star-forming phase 
in spheroids is completed at z > 1. 

Note that, since the clustering amplitude strongly depends on the quantity Mhaio/Msph; measurements 
of the angular correlation function w(9) are in principle able to discriminate amongst different models 
of SCUBA galaxies and in particular to determine both their star-formation rate, via the amount of 
baryonic mass actively partaking the process of star formation, and the duration of the star-formation 
burst. 

A first attempt to measure the angular correlation function of S > 5 mJy SCUBA sources has been 
recently presented by [O]. Although such measurements are dominated by noise due to small-number 
statistics, it is nevertheless interesting to note that - as illustrated by the left-hand panel of Figure 2 
- our model (with a preference for the -/Wh a io/-^sph = 100 case) shows full consistency with the data 
(kindly provided by S. Scott). 

Another possible way to obtain some information on the nature of SCUBA sources via their 
clustering properties is provided by the predictions of p| for Lyman-break galaxies to be the low- 
to intermediate-mass tail of primeval spheroidal galaxies, with Tg ~ 1 — 2 Gyr and a star-formation 
rate ranging from a few to a hundred M@ yr _1 . In Figure 2 (right-hand panel) we then plotted the 
predicted w{6) for those sources with S > 1 mJy (corresponding to M sp h ^ 1O 1O M ), expected to be 
found within the redshift range 2.5 < z < 3.5, covered by the original (Steidel et al., 1996) sample. As 
in the former case, the higher curve is for Mh a i /M sp h = 100, the lower one for Mhaio/Msph = 10. The 
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Figure 2: Left-hand panel: angular correlation function of S > 5 mJy SCUBA sources. Model 
predictions are illustrated by the dashed lines (where the higher one is for Mhaio/Mgph = 100 and 
the lower one for M^io/M^^ = 10), while data points represent the |l3| measurements. Right-hand 
panel: angular correlation function of LBGs. Model predictions are shown by the dotted lines (the 
higher one for Mh a i /M sp h = 100, the lower one for Mh a i /M S ph = 10), while data points represent the 
[Si measurements. 



data points show the Q measurements. Even though large errors once again affect the observational 
findings, it is nevertheless clear that the predicted trend for w(9) can correctly reproduce the data 
for high (80 — 100) values of the Mh a i /M sp h ratio. This result is consistent with the predictions by 
[f§ and implies a well defined relationship between SCUBA galaxies and LBGs. Furthermore, it also 
confirms the expectations for a small fraction (on the order of a few percent) of the total mass to be 
confined into stars. 

Finally, it is also worth noticing that our predictions are in agreement with the strong clustering 
of EROs recently detected by [||, since we expect ([||) these objects to be the direct descendants of 
SCUBA galaxies, and therefore to exhibit the same clustering properties. 



2 Power Spectrum of Temperature Fluctuations 

An issue intimately connected with the analysis of galaxy clustering is the study of the contribution 
of unresolved sources (i.e. sources with fluxes fainter than some detection limit Sd) to the background 
intensity. Its general expression is given by: 
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(see e.g. ||]), where dN/dS denotes the differential number counts, L max and L min are respectively 
the maximum and minimum local luminosity of the sources, K(L, z) is the K-correction, z max is the 
redshift when the sources begin to shine, z(Sd, L) is the redshift at which a source of luminosity L is 
seen with a flux equal to the detection limit Sd, $(L, z) is the luminosity function (i.e. the comoving 
number density of sources per unit (ilogL), and x is the comoving radial coordinate. 

The intensity fluctuation 51 due to inhomogeneities in the space distribution of unresolved sources 
is then given by eq. (|J), with the quantity $>(L,z) replaced by #$(L, z). It is easily shown that 
the angular correlation of such intensity fluctuations C{9) = (61(9', (f)') 51(9" , eft")), where (9',<f>') and 
(9" ', 4>") define two positions on the sky separated by an angle 9, can be expressed as the sum of two 
terms Cp and Cc, the first one due to Poisson noise (i.e. fluctuations given by randomly distributed 
objects), and the second one owing to source clustering. It is possible to show (J7j) that, in the case of 
highly clustered sources, the Poissonian term Cp is negligible with respect to the one due to clustering. 
In the following we therefore only concentrate on temperature fluctuations caused by the Cc term 
(hereafter simply called C). 

By making use of the quantities previously defined and of eq. (y), the clustering term C takes the 
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as) = (s) C . , di 6 - <M * z) irw (I) f *• «■-• «>• (4) 

where the effective volume emissivity j e g is expressed as: 

/■min[L maXiL(s } ] 

JoS= <!>(L,z)K(L,z)LdlogL (5) 

C(0) in eq.(Q) has been evaluated separately for the three cases of low-, intermediate- and high- 
mass objects, by plugging in eq. (|5|) the appropriate expressions for the luminosity function. The total 
contribution of clustering to intensity fluctuations has then been derived by adding up all the values 
of C{6) obtained for the different mass intervals and by also taking into account the cross-correlation 
terms between objects of different masses, according to the expression 

c TOT {e) = £ y/cmctf), (6) 

where the indexes i,j stand for high, intermediate and low masses. 

Note that, the quantity 6g ff (M m i n , z) in eq.@ should indeed be read as 6g ff (M m i n , M max , z), where 
M max is the maximum halo mass corresponding to the maximum visible bulge mass (i.e. upper limit 
for the mass locked into stars), which corresponds in eq. (J2|) to a replacement in the upper limit of the 
integrals of oo with Af max . 

The angular power spectrum of the intensity fluctuations can then be obtained via (see Q): 

ST(9) = ((AT) ) = — ^exp ^—j - ij exp (- — j / (^— j , (7) 

which relates intensity and temperature fluctuations. Figure || shows the predicted values for the 
quantity STi = yjl(l + l)Ci/2ir (in units of K) at respectively 353 GHz (850 /xm - left-hand panel) 
and 545 GHz (550 ^m - right-hand panel) - the central frequencies of two of the channels of the 
High Frequency Instrument (HFI) of the ESA's Planck mission - as a function of the multipole I up 
to I = 1000. Results are plotted for two different values of the source detection limit (Sd = 100 and 
10 mJy for the 353 GHz case and Sd = 450 and 45 mJy for the 545 GHz case) and the usual two values 
of Mhaio/Msph- Also shown, for comparison, is the power spectrum of primary (CMB) anisotropies 
(solid line) predicted for the cosmology specified in the caption, computed with the CMBFAST code 
developed by |T^| . 

At 850 /mi, our model predicts fluctuations of amplitude due to clustering comparable to (and 
possibly even larger than) those obtained for primary CMB anisotropies at I ^ 50. This is because most 
of the clustering signal comes from massive galaxies with fluxes S ^ 10 mJy, which lie at substantial 
redshifts and are therefore highly biased tracers of the underlying mass distribution. Also the strongly 
negative K-correction increases their contribution to the effective volume emissivity [eq. (||)] and 
therefore to the fluctuations. The 550 /um case is even more striking, since the contribution from 
clustering is expected to be more than an order of magnitude greater than the one originating from 
primordial fluctuations, regardless of the flux detection limit. 

This implies that important information on the clustering properties of faint sub-mm galaxies (and 
hence on physical properties such as their mass and/or the amount of baryons involved in the star- 
formation process) will reside in the Planck maps at frequencies greater than 353 GHz where, however, 
the dominant signal is expected to come from interstellar dust emission. In order to show this effect, 
in Figure 3 we have also plotted the expected contribution from galactic dust emission averaged all 
over the sky (upper dashed-dotted curves). This was derived from the IRAS maps at 60 /im, rescaled 
at the frequencies under exam by assuming a grey-body spectrum with Tdjjst = 18 K. As already 
anticipated, this signal appears to be the dominant one in both the 353 and 545 GHz Planck channels. 
Nevertheless, it is still possible to extract information on the nature of sub-mm galaxies if one restricts 
the analysis to high galactic latitude regions (i.e. 80° < b < 90°, lower dashed-dotted curves in Figure 
3), which are the least affected by galactic dust emission. In fact, as it can be seen from Figure 3, 
the dust contribution in this region becomes less important than the one due to the clustering of 
unresolved sources for Z ^ 80. 
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Figure 3: Predicted power spectrum of temperature fluctuations 8Ti = y/l(l + l)Ci/2ir (in units of 
K) as a function of the multipole I at 353 GHz (left-hand panel) and 545 GHz (right-hand panel), 
the central frequencies of two Planck/HFI channels. Left-hand panel: dashed lines are for a detection 
limit Sd = 100 mJy, dotted ones for Sd = 10 mJy. In both cases higher curves are obtained for 
Mhaio/M sp h = 100, lower ones for Mh a ioA4ph=10- The solid line represents the power spectrum of 
primary CMB anisotropics as predicted by a standard Cold Dark Matter model for a ACDM cosmology 
(A = 0.7, Qq = 0.3, /io = 0.7). Right-hand panel: as in the former case but with flux detection limits 
of 450 (dashed lines) and 45 (dotted lines) mJy. The overlap observed here is due to paucity of bright 
(S ;> 50 mJy) sources. In both panels the upper dashed-dotted curve represents the contribution from 
galactic dust emission averaged all over the sky, while the lower one illustrates the case obtained by 
restricting the analysis to high galactic latitudes (80° < b < 90°). 
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